Electronic and spin transport properties of zigzag graphene nanoribbon mediated by 

metal adatom: QUAMBO-NEGF approach 
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Electronic and spin transport properties of zigzag graphene nanoribbon (ZGNR) with alkali-, 
group-Ill and 3d-transition-metal adatoms on the hollow center are investigated by quasi-atomic 
minimal basis orbits (QUAMBOs), a first-principles tight binding (TB) scheme, combined with non- 
equilibrium Green's function (NEGF). QUAMBOs-TB parameters of carbon atoms in the middle 
scattering region are constructed by divide-and-conquer approach. It is found that the transmission 
is generally reduced by metal adatom, which can be well explained by the geometry of 7r-orbital 
and the corresponding interaction strength between 7r-orbital and metal adatom. Naturally, the 
transport at the Fermi level is not sensitive to metal adatom because of weak interaction between 
the localized edge state and metal adatom. Spin transport near the Fermi level is induced by Co, 
Fe and V adatom, and higher than 50% transmission spin polarization with a board energy window 
around 0.34eV, 0.36eV and 0.56eV respectively. Thereinto, a higher transmission of one spin state 
is originated from a less distortion of the corresponding 7r-orbital by metal adatom rather than high 
DOS of this spin state. Further the effect of the width of ZGNR on electronic and spin transport 
properties is discussed. 

PACS numbers: 72.80.Vp, 73.20.At, 31.15.A-, 
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I. I. INTRODUCTION 

Graphene [l| , a true two-dimensional crystalline lattice 
formed by carbon atoms, has attracted tremendous inter- 
est in physics, chemistry, and material science. The novel 
electronic structures of graphene with a linear dispersion 
relation near the Dirac point generate unusual transport 
properties [l|-0] , such as the conductivity linear with the 
gate voltage resulting from the linear spectrum 0, Q , un- 
conventional integer quantum Hall effect [H, Hj] , and high 
carrier mobility |5J which make graphene as a good can- 
didate in electronic devices. In addition, there have been 
a lot of interests in modifying the electronic, spintron- 
ics, transport, as well as magnetic properties of graphene 
through adsorption of various atoms and molecules on 
graphene interacting with graphene sheet (^-tllj. Simi- 
lar to the proposal of biosensor by Schedin et al. [i~2j . 
it was observed gate-tunable spin transport in graphene 
at room temperature fl3l - FT5j . First principle investiga- 
tion of the effects of metal contact on graphene [H, [I?} 
was stimulated by the metallic contacts in the transport 
measurement. On the other hand, the effect of adatom, 
dimmer and atom-chain of transition metal on electronic 
and spin transport had also been explored [l8l - [20j . 

Due to the huge computational demand for first prin- 
ciple calculation of large scale system, lot of efforts had 
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been made to extrac t tig ht binding parameters from first 
principle calculation |l7l.l2l| . Recently quasi-atomic min- 
imal basis set orbitals (QUAMBOs) was developed by 
Yao et al., in which TB Hamiltonian and overlap matri- 
ces of a given system can be accurately extracted from ab 
initio calculations without any fitting procedures [22l.[23|. 
Furthermore, QUAMBOs scheme enables the applica- 
tion of divide-and-conquer approach benefitted from the 
dependence of local environment of QUAMBOs-TB pa- 
rameters as shown in Rcfs. [IH H3]- On this basis, 
non-equilibrium Green's function (NEGF) method com- 
bined with Landauer' formalism |25l - [27| are preformed 
to calculate the electronic transport. Therefore such a 
QUAMBO-NEGF scheme can be applied to study the 
electronic transport in l arg e atomistic systems without 
sacrificing the accuracy |24j . 

Our previous work investigated the bond and charge 
transfer between various metal adatom, including alkali-, 
group-Ill and 3d-transition-metal adatom, and graphene 
system [l(| • In this letter, we further investigate the 
electronic and spin transport through ZGNR in the pres- 
ence of these metal adaoms by QUAMBO-NEGF scheme. 
There are two distinct differences between our work and 
previous similar works [l8l . |20j . One is that spin un- 
polarized basis is adopted in our first-principles calcu- 
lation of ideal ZGNR [24], HH, since the ground state 
of ZGNR under spin polarized basis is not stable at fi- 
nite temperature. Certainly, for graphene with metal 
adatom as shown in Fig. 1, spin polarized basis should 
be used in first-principles calculation The other 



is that the Fermi level of the middle scattering region 
is shifted to zero [HI, |2(J. However, here TB param- 
eters near the interfaces equal to the counterparts at 
the leads, in order to be consistent with the screen- 
ing of metal adatom at the interface, and therefore the 
Fermi level of the middle scattering region may not be 
shifted to zero. Though above differences and different 
width of ZGNR, the trends of electron current with a 
certain dominant spin state [2(| are consistent with 
our work. We find that the current with a certain dom- 
inant spin state is originated from a less distortion of 
the corresponding 7r-orbit by metal adatom. Furhtcr- 
more, a divide-and-conquer approach is adopted to con- 
struct QUAMBOs-TB parameters for ZGNR with metal 
adatom on hollow center, which goes beyond the con- 
struction of QUAMBOs-TB parameters in the presence 
of one type of atom [28[ . 

This paper is arranged as follows. In Sec. 2, we present 
the VASP calculation details and formula of the trans- 
mission and density of states based on QUAMBO-TB 
Hamiltonian and overlap matrices; in Sec. 3, we describe 
the reconstruction of TB parameters for zigzag graphene 
nanoribbon (ZGNR) with various metal adatoms and in- 
vestigate the effect of metal adatom and the width of 
ZGNR on electronic and spin transport properties; and 
finally we give a brief summary in Sec. 4. In our investi- 
gation, the transmission is reduced by metal adatom and 
is well explained by density of states, the distortion of 
QUAMBOs-TB parameters by metal adatom, the geom- 
etry of 7r-orbit and then the interaction between 7r-orbit 
and metal adatom. in sections V and VI. 



II. II. COMPUTATIONAL DETAILS 

The first-principles calculations are performed based 
on the density functional theory (DFT) with gener- 
alized gradient approximation (GGA) in the form of 
PBE (Perdew-Burke-Ernzerhof) [29j implemented in the 
VASP [MS] (Vienna Ab initio Simulation Package) 
code, including spin polarization and dipole moment cor- 
rections [33l . |34 |. Valence electrons are treated explic- 
itly and their interactions with ionic cores are described 
by PAW (Projector Augmented Wave) pseudopotentials 
[35l l36j . The wave functions are expanded in a plane wave 
basis set with an energy cutoff of 600 eV. A k-point sam- 
pling of 6x6x1 Monkhorst-Pack grids in the first Bril- 
louin zone of the supercell and a Gaussian smearing with 
a width of 0.05 eV is used in the calculations. All atoms 
in the supercell are allowed to relax until the forces on 
each atom are smaller than 0.01 eV/A. The supercell 
dimensions are kept fixed during the relaxation. 

The adatom/graphene system is modeled by having 
one adatom in a 4x4 parallelogram graphene supercell 
with periodic boundary conditions as shown in Fig. 1. 
The primitive cell of graphene is a parallelogram with 
two carbon atoms. The lattice constant obtained from 
our calculation is 2.46A, agrees well with experimental 
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FIG. 1: (Color online) Three adsorption sites for an adatom 
on graphene: hollow center (H), bridge center (B), and top 
(T) sites in a periodic 4x4 graphene supercell. In first- 
principles calculation, H site is stable for metal adatom ad- 
sorbed on the graphene. 



value. The dimension of the supercell in the z direction 
is 15A, which allows a vacuum region of about 12A to 
separate the atoms in the supercell and their replicas. 
Usually there are three typical sites for adatoms posi- 
tioned on graphene at the top of a carbon atom, labeled 
top (T) site, at the middle of a carbon-carbon bond, la- 
beled bridge (B) site, and at the hexagonal center site, 
labeled hollow (H) site, respectively, as indicated in Fig. 
1. In our calculations, H site is stable for metal adatoms 
adsorbed on the graphene under study. 

Using the self-consistent wave functions from the 
plane-wave basis VASP calculations, a set of structure- 
dependent quasi-atomic minimal basis orbitals (QUAM- 
BOs) for each adatom/graphene system is constructed 
following the scheme developed previously [22j,|23j. While 
highly localized on atoms and exhibiting shapes close to 
the orbitals of the isolated atoms, the QUAMBOs span 
exactly the same occupied subspace as the wave func- 
tions determined by the first-principles calculations with 
a large plane- wave basis set. The Hamiltonian and over- 
lap matrices in the QUAMBOs representation give ex- 
actly the same energy bands and wave functions of the 
occupied electronic states as those obtained by the fully 
con verged first-principles calculations using the VASP 
[HI [22, 23]. Therefore, using the wave functions from the 
VASP calculations, accurate tight-binding (TB) Hamil- 
tonian and overlap matrix elements are calculated by 
QUAMBO approach [H |23|. The electronic transmis- 
sion probability in such systems can be calculated as: 

T a (E) = Tr[T L AE)G*{E)T R , a {E)G*{E)], (1) 

R A 

T L>(T (E) = i[J2(E)-J2( E )}- (2) 

L.cr L,(T 

Here Tl i(J (E) is the coupling function of left semi-infinite 
lead (L, ZGNR on the left of Fig. 2b), and J2t.A E ) 
is the conjugate of the advanced self energy Y^Loi-®). 
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TliL a(E) is calculated through surface Green's function 
iteratively pj - E^ . 

R 

^2(E) = (H LM a - ES LM a)9L,a ( H LMc7 _ ES LMa) ' 

(3) 

In the present calculation, whether the spin state indexed 
by a is degenerate or not depends on the adsorbed metal 
adatom on ZGNR. M and R stand for the middle scat- 
tering region and the right lead respectively, a is the 
surface Green's function of the left lead. Hlmo {Slmu) 
is the Hamiltonian (overlap) matrix between left lead and 
the middle scattering region. The same calculation can 
also be performed for the right lead to obtain YIr a {E)- 
The retarded Green's function is calculated by 

(R R \ 1 

Es M ,„ - h m><7 -J2( E n > 
L,a R,a J 

(4) 

where Hm,u [Sm,<t) is the Hamiltonian (overlap) matrix 
of the middle scattering region. E is the energy for elec- 
tron injecting from left reservoir. E = implies that the 
Fermi energies of leads have been shifted to zero, while 
the Fermi energy of the middle scattering region has not 
been shifted to zero due to the presence of metal adatom. 
We stress that the conductance of system is under zero 
bias between left and right leads. Meanwhile, the den- 
sity of states (DOS) p a (E) = -ImTr [G^{E)S M ,a] /V in 
open boundary system of lead-middle scattering region- 
lead should be different from DOS in periodic boundary 
system by first-principle calculation. 

III. III. RESULT AND DISCUSSION 

A. Reconstruction of TB parameters for ZGNR 
with metal adatom on the hollow center by 
divide-and-conquer approach 



by metal adatom, except for those of six carbon atoms 
close to metal adatom, as shown in Fig. 3(b). The local 
environment dependence enables us to easily construct 
a wider system and estimate the effect of the width of 
ZGNR lead on transport property. 
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FIG. 2: (Color online) (a) The periodic 4x4 graphene su- 
percell shown in Fig. 1 is arranged into a rectangle shape, 
(b-c) The middle scattering region in the blue frame, i.e., the 
conductor part, is composed by zigzag graphene nanoribbon 
(ZGNR) with an metal adatom at the hollow center. Passi- 
vating hydrogen atoms at the edges of ZGNR are not shown. 
The primitive unit (PU) of ZGNR lead is in the red frame. 
The width of ZGNR W is 8 and 12 in (b) and (c) respectively. 

Magnetism in graphene is induced by 3d transition 
metal, such as Co, Fe and V. In this case, TB parameters 
indexed by spin up and down indices are not degener- 
ate due to the induced magnetism in graphene sheet, as 
an example in the presence of Co adatom shown in Fig. 
3(b). Thereinto, spin splitting effect is strongest for six 
carbon atoms closest to these metal adatoms as shown 
in Fig. 3(a). However, this kind of spin splitting ef- 
fect is slight for those carbon atoms far away from these 
metal adatoms. Therefore it is reliable to retain spin- 
degenerate QUAMBO-TB parameters for those carbon 
atoms close to the zigzag edge in the middle scattering 
region. 



In order to study the transport property of zigzag 
graphene nanoribbon (ZGNR) with various metal 
adatom adsorbed at the hollow center, we choose a gen- 
eral structure of lead-middle scattering region-lead as 
shown in Fig. 2(b,c). The lead and the middle scattering 
region are a semi-infinite ZGNR and a ZGNR with metal 
adatom respectively. To construct the middle part, we at 
first reshape the structure as shown in Fig. 1, into a rect- 
angle as shown in Fig. 2(a), and then replace the center 
part of ZGNR in the middle scattering region by this 
rectangle. The feasibility of this reconstruction scheme 
is assured by the local environment dependence of TB pa- 
rameters obtained by QUAMBOs as the example of car- 
bon atom chain between ZGNR was represented in Ref. 
[28| . Here we check the screening effect of metal adatom 
at lead-middle scattering region interface and find that 
the TB parameters of carbon atoms are slightly distorted 



B. Electronic and spin transport properties in the 
presence of metal adatom and comparison with 
some other published works 

For ideal ZGNR with W=8 (here W being the number 
of carbon atoms in one arrow as shown in Fig. 2b), the 
transmission is quantized, as shown in Fig. 4(a). How- 
ever the transmission around a narrow region of Fermi 
level with the width of 0.04eV is larger than 1, similar 
to the transmission in Ref. [13] ■ QUAMBOs-TB parame- 
ters are reliable as one perfect transmission for the energy 
range [-1.06eV, 1.24eV] with a width of 2.30eV is well 
reproduced, as shown in Fig. 4(a). We plot the spatial 
density of states (SDOS) of each carbon atom at E=-l, 
-0.6, 0, 0.6 and 1 eV in Fig. 5(a-e). It is found that high 
SDOS horizontally distributes at the center of ZGNR as 
E is less than 0, and high SDOS moves toward the zigzag 
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edge little by little and finally an edge state occurs when 
E approaches 0. On the contrary, SDOS vertically dis- 
tributes when E is larger than 0. Therefore the geometry 
of 7r-orbit determines the interaction strength between 
7r-orbit and metal adatom and then transport properties 
through ZGNR. 
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FIG. 3: (Color online) (a) Ordering for carbon atoms in rect- 
angle shape as shown in Fig. 2(a). (b) On-site energies of 
carbon atom in the presence of Co and K adatom. 7r-orbits 
of six carbon atoms closest to the metal adatom (i.e. car- 
bon # 10, 11, 14, 15, 18 and 19) are strongly affected by 
the metal adatom, while 7r-orbits of other carbon atoms are 
slightly distorted. The local-environment dependence of TB 
parameters enables the reconstruction of middle scattering 
region as shown in Fig. 2. 

The transport properties of ZGNR with various metal 
adatom are shown in Fig. 4(a-b,d-f). It is found that the 
transmission may be spin-independent as shown in Fig. 
4(a-b) and spin-dependent as shown in Fig. 4(d-f), which 
depends on the net magnetism induced by metal adatom. 
First of all, we discuss the spin- independent transport. It 
is found that the transmission within the energy range [- 
2.0eV, -0.02eV] is greatly reduced by the distortion of 
metal adatom. Thereinto, the perfect channel within [- 
1.06eV, -0.02eV] for ideal ZGNR with W=8 is blocked, 
which is embodied in localized SDOS at the conductor- 
lead interfaces as shown in Fig. 5(f-g), as 7r-orbit with 
high SDOS distributing horizontally in the center inter- 
acts strongly with metal adatom on hollow center. How- 
ever, the perfect channel within [0.02eV, 1.24eV] remains 
a finite transmission in the presence of metal adatom, 
since low SDOS is in the vicinity of metal adatom and 
then a weak interaction occurs between metal adatom 
and 7r-orbit as shown in Fig. 5(i-j). It can be inferred 
that the finite transmission in this energy range is orig- 
inated from the unchanged geometry of SDOS, though 



SDOS is distorted by metal adatom. It is well under- 
stood that the transmission at the Fermi level E=0 is not 
altered in the presence of various metal adatom, since the 
edge state close to zigzag edge is far away from the metal 
adatom and then weakly interacts with metal adatom. 
Its corresponding SDOS shows that the edge state has 
not been distorted by the metal adatom as shown in Fig. 
5(h). 

It is interesting that the transmission curves collapse 
as the energy ranges from [-2.0eV, 1.24eV] in the presence 
of Al, In, K and Na, as shown in Fig. 4(a). The trans- 
mission curve within the energy range [-2.0eV, 1.24eV] 
in the presence of Ni adatom almost collapses with those 
mediated by Al and K adatoms, as shown in Fig. 4(b). 
Different from the above four metal adatoms, Ni adatom 
provides two transmission peaks at -0.94eV and -0.34eV 
as a result of peaks in high density of states (DOS), as 
shown in Fig. 4(c). As the energy is lower than -1.75 
eV, the transmission is not sensitive to metal adatoms, 
since 7r-orbit has a weak interaction with metal adatom. 
As the energy is higher than 1.24eV, spin-independent 
transport mediated by metal adatom shows three differ- 
ent trends as shown in Fig. 4(a) and 4(b). This is induced 
by the types of metal adatoms, (1) akali-metal Al and 
In, (2) group-Ill metal K and Na, and (3) 3d transition- 
metal Ni. Their corresponding outer electron structures 
are 3s 2 3p 1 , 5s 2 5p 1 , 3s 1 , 4s 1 and 3d 8 4s 2 respectively. It 
shows that the transport at high energy is determined by 
the outer electron structure of metal adatoms. We should 
stress that Ni/graphene system has no net magnetic mo- 
ment due to charge transfer from Ni adatom to graphene, 
though isolated Ni atom has of magnetic moment fioj . 
Consequently, the transmission is spin-degenerate in the 
presence of Ni adatom on ZGNR. 

The spin splitting in transport is obvious within the 
energy range near Fermi level once there is magnetism 
induced in graphene sheet by transition metal adatom, 
i.e. Co, Fe and V. However, no spin splitting presents in 
transmission for the energy far away from the Fermi level, 
which applies that magnetism mainly originates from n- 
orbit consistent with DOS from first-principles calcula- 
tion. In the presence of Co adatom, the transmission 
for electron with spin down (Tdown) within the energy 
range [0.02, 1.0] eV is higher than the transmission for 
electron with spin up (Tup). Meanwhile, its density of 
states (DOS) for electron with spin down is lower than 
the opposite within the above energy range, as shown 
in Fig. 4(d). Therefore, a high transmission for electron 
with spin down is due to less distortion of 7r-orbit indexed 
by spin down. This conclusion is also embodied in the 
less distortion of the on-site energy indexed by spin down 
in the presence of Co adatom as shown in Fig. 3(b). The 
relation between the distortion and the transmission in 
the above applies to all non-degenerate spin transport in 
the presence of other transition metal adatoms Co and V 
as shown in Fig. 4(e-f). It is found that the transmission 
curve within the energy range [-2, 1.24] eV collapses with 
that mediated by Al adatom, when the spin state results 



FIG. 4: (Color online) Transmission versus the energy 
through ZGNR in the presence of (a) Al, In, K and Na adatom 
and (b) Al, K and Ni adatoms adsorbed on the hollow center 
of ZGNR. Curves of the transmission versus the energy for 
exact versus reconstructed ZGNR are plotted in (a). "RC, 
W8-lead" stands for reconstructed ZGNR HJHi]. ( c ) Den- 
sity of states (DOS) versus the energy for ZGNR with Al, In, 
K, Na and Ni adatoms. (d-f) show transmission for electron 
with spin-up and spin-down indices and the corresponding 
DOS versus the energy through ZGNR in the presence of Co, 
Fe and V adatom respectively. Transmission spin polariza- 
tion (TSP) in the presence of Co, Fe and V adatom is shown 
in the right inset of (d). The transmission in the presence of 
Al adatom is plotted in (d-f) for comparison. The width of 
ZGNR is W=8. 
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FIG. 5: Spatial density of states (SDOS) of each carbon atom 
in (a-e) pure ZGNR and (f-j) ZGNR with Al metal adatom 
on hollow center. SDOS mainly comes from 7r-orbit of carbon 
atoms. The width of ZGNR W is 8. 



could not induce a finite transmission around the Fermi 
level [2(J. Though different basis adopted in calculation 
and different width of ZGNR from our work was cho- 
sen, it was found that the electron with spin down (spin 
up) in the presence of Co (Fe) adatom provides a higher 
transmission within some certain energy ranges [Tel l20j|. 
which is consistent with our work. 



in a less distortion to 7r-orbit, as shown in Fig. 4(d-f). We 
also investigate the transmission spin polarization (TSP), 
which is defined as (Tup-Tdown) / (Tup+Tdown) . Cur- 
rent with spin up (down) dominates when TSP is pos- 
itive (negative). In view of the transmission and TSP, 
the energy range with a range width of 0.34eV, 0.36eV 
and 0.56eV above the Fermi level is suitable for spin ma- 
nipulation mediated by Co, Fe and V respectively, where 
the absolute value of TSP is larger than 0.5. As we dis- 
cuss above, the transmission at the energy higher than 
1.24eV sensitively depends on the outer electron struc- 
ture of metal adatom. 

As we mentioned in the introduction, in other similar 
works on transport property of ZGNR in the presence of 
metal adatom, the Fermi level of the middle scattering re- 
gion is shifted to zero [H, Ho| . Furthermore, ZGNR lead 
has a gap when spin polarized basis is chosen, therefore 
the transmission near the Fermi level is zero for ZGNR 
with Co adsorbed on the hollow center [Tsj . Similarly, 
it was found that a single Fe atom adsorbed on ZGNR 



C. Effect of the width of ZGNR lead on electronic 
and spin transport 

The trend of the effect of ZGNR lead on electronic 
and spin transport is highly desired, as electronic devices 
are usually wider than the system under study. How- 
ever, it is difficult to include a wide ZGNR lead in first- 
principles calculation due to time and memory consum- 
ing. We continue to adopt divide- and- conquer approach 
to construct a ZGNR lead with W=12. The comparison 
of the transmission for two pure ZGNR leads is shown 
in Fig. 6(a). The transmission is quantized for pure 
leads, where the quantity of transmission depends on the 
number of eigenchannels available for the certain energy. 
Since the wider lead provides more channels for electronic 
transport, therefore the energy range for each quantized 
transmission becomes narrower as the width of ZGNR 
increases. 

ZGNR with Co adatom on the hollow center is chosen 
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FIG. 6: ((Color online) (a) Transmission versus energy 
through ZGNR leads with W=8 and W=12. (b) The trans- 
mission and transmission spin polarization (TSP, shown in 
the inset) through ZGNR in the presence of Co adatom as 
the width of ZGNR W is 8 and 12 respectively. 



as an example to estimate the effect of ZGNR on elec- 
tronic and spin transport as shown in Fig. 6(b). The 
curve of the transmission versus the energy in the case 
of W=12 ZGNR lead is similar to the counterpart in the 
case of W=8 ZGNR lead, in which the transmission in- 
creases with the absolute value of the energy except for 
those energy ranges with no transmission at all. The dif- 
ference between two cases is that wider ZGNR provides 
more channels for electronic transport and therefore the 
transmission is generally higher, except for those energy 
ranges with no transmission at all. We also investigate 
the behavior of TSP versus the energy and find that the 
energy range corresponding to a finite TSP becomes nar- 
rower and closer to the Fermi level as the width of ZGNR 
increases. Similarly, taken into account a finite trans- 
mission, the energy range [0.3, 0.52] eV is suitable for 



spin manipulation, in which the absolute value of TSP 
is larger than 0.5, as the width of ZGNR lead is W=12. 
It is found that not only the effect of spin splitting but 
the effect of width of ZGNR lead on the transmission 
at E—Q is negligible. As we discussed above, this phe- 
nomenon occurs because the edge state localized at the 
zigzag edge is slightly distorted by the metal adatom. 



IV. IV. SUMMARY 

Electronic and spin transport properties of metal 
adatom adsorbed on zigzag graphene nanoribbon 
(ZGNR) are investigated by QUAMBO-NEGF combined 
with divide-and-conquer approach. The effect of metal 
adatom on the transmission has been investigated. Non- 
degenerate spin transport near the Fermi level is induced 
by transition metal adatom such as Co, Fe and V. Spin 
manipulation has also been analyzed in term of transmis- 
sion spin polarization. Further the trend of the electronic 
and spin transport properties through a wider ZGNR 
lead is discussed. 
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